Linear-response theory of spin Seebeck effect in ferromagnetic insulators 
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We formulate a linear response theory of the spin Seebeck effect, i.e., a spin voltage generation from 
heat current flowing in a ferromagnet. Our approach focuses on the collective magnetic excitation 
of spins, i.e., magnons. We show that the linear-response formulation provides us with a qualitative 
as well as quantitative understanding of the spin Seebeck effect observed in a prototypical magnet, 
yttrium iron garnet. 



I. INTRODUCTION 

The generation of spin voltage, i.e., the potential for an 
electron's spin to drive spin currents, by a temperature 
gradient in a ferromagnet is referred to as the spin See- 
beck effect (SSE). Since the first observation of the SSE 
in a ferromagnetic metal NisiFeigjfi this phenomenon has 
attracted much attention as a new method of generat- 
ing spin currents from heat energy and opened a new 
possibility of spintronics devices^ The SSE triggered the 
emergence of the new field dubbed "spin caloritronics"^ 
in the rapidly growing spintronics community. Moreover, 
as the induced spin voltage can be converted into electric 
voltage through the inverse spin Hall effect 5 at the at- 
tached nonmagnetic metal, this phenomenon put a new 
twist on the long and well-studied history of thermoelec- 
tric research^ 

One of the canonical frameworks to describe nonequi- 
librium transport phenomena is linear-response theory J 
Having been applied to a number of transport phenom- 
ena, linear-response theory has been so successful be- 
cause it is intimately related to the universal fluctuation- 
dissipation theorem. Up to now, however, the linear- 
response formulation of the SSE has not been known 
mainly because, unlike the charge current, the spin cur- 
rent is not a conserved quantity. Therefore, it is of 
great importance to formulate the SSE in terms of linear- 
response theory. 

Concerning the SSE, a big mystery is now being estab- 
lished, which is, how can conduction electrons sustain the 
spin voltage over such a long range of several millimeters 1 
in spite of the conduction electrons' short spin-flip diffu- 
sion length, which is typically of several tens of nanome- 
ters? A key to resolve this puzzle was reported by a re- 
cent experiment on electric signal transmission through 
a ferromagnetic insulator 8 which demonstrates that the 
spin current can be carried by the low-lying magnetic 
excitation of localized spins, i.e., the magnon excitations, 
and that it can transmit the spin current as far as several 
millimeters. Subsequently, the SSE was reported to be 
observed in the magnetic insulator LaY2Fe50i2 despite 
the absence of conduction electrons.— These experiments 
suggest that contrary to the conventional wisdom over 
the last two decades that the spin current is carried by 



conduction electrons^ the magnon is a promising can- 
didate as a carrier for the SSE. 

The purpose of this paper is twofold. First, we analyze 
the SSE observed in LaY2FesOi2^ (hereafter referred to 
as YIG) in terms of magnon spin current, i.e., a spin cur- 
rent carried by magnon excitations. Second, we develop 
a framework for analyzing the SSE by means of the stan- 
dard linear-response formalism 7 which is amenable to the 
language of the magnetism community M This allows us 
to describe the spin transport phenomena systematically, 
and it can be easily generalized to a situation including 
degrees of freedom other than magnons, e.g., conduction 
electrons and phonons, to describe a more complicated 
process in the case of metallic systems^ 

The plan of this paper is as follows. In Sec. |Hl we 
present a linear-response approach to the "local" spin in- 
jection by thermal magnons, in which the spin injection 
is driven by the temperature difference between the fer- 
romagnet and the attached nonmagnetic metal. Next, in 
Sec. IIIII we develop a linear-response theory of the "non- 
local" spin injection by thermal magnons, in which the 
spin injection is driven by the temperature gradient in- 
side the ferromagnet. As one can see below, this process 
can explain the SSE observed in YIGi^ Finally, in Sec. lIVI 
we summarize and discuss our results. 



II. "LOCAL" SPIN INJECTION BY THERMAL 
MAGNONS 



We start by briefly reviewing the SSE experiment for 
YIG. 9 Figure [T] shows the experimental setup where sev- 
eral Pt terminals are attached on top of a YIG film in 
a static magnetic field H$z (3> anisotropy field) which 
aligns the localized magnetic moment along z. A tem- 
perature gradient VT is applied along the z axis, and it 
induces a spin voltage across the YIG/Pt interface. This 
spin voltage then injects a spin current / s into the Pt ter- 
minal (or ejects it from the Pt terminal). A part of the 
injected/ejected spin current / s is converted into a charge 
voltage through the so-called inverse spin Hall effect £ 
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FIG. 1: Experimental setup for observing the SSE.— Inset: 
Schematics of the spatial profile of the observed voltage. 

where |e|, 0#, p, and w are the absolute value of electron 
charge, spin Hall angle, resistivity, and width of the Pt 
terminal (see Fig. [T]), respectively. Hence, the observed 
charge voltage Vishe is a measure of the injected/ejected 
spin current I s . 

To investigate the SSE observed in YIG, we consider 
a model shown in Fig.[2ja). While YIG is a ferrimagnet, 
we model it as a ferromagnet since we are interested in 
the low-energy properties. The key point in our model 
is that the temperature gradient is applied over the in- 
sulating ferromagnet, but there is locally no temperature 
difference between the ferromagnet and the attached non- 
magnetic metals, i.e., Tn 1 = Tp 1 = Ti, T^ 2 = Tp 2 = T2, 
and Tn 3 = Tp 3 = T3. We assume that each domain 
is initially in thermal equilibrium without interactions 
with the neighboring domains, and then calculate the 
nonequilibrium dynamics after we switch on the interac- 
tions. Note that this procedure is essentially equivalent 
to that used by Luttinger— to realize the initial condition 
mentioned above. 

Let us consider first the low-energy excitations in the 
ferromagnet. In the following, we focus on the spin- 
wave region where the magnetization M(r) fluctuates 
only weakly around the ground state value M s z with 
the saturation magnetization M g , and we set M/M s = 
(1 — m 2 /2)z + m to separate the small fluctuation 
part m (_L z) from the ground state value. Then, 
the low-energy excitations of M are described by bo- 
son (magnon) operators a q and a q through the rela- 
tions^ m q = y/TJSoa_ q and m q = y/TJ Soa q where 
m = (m x ± im v )/^/2, Sq is the size of localized spins, 
and m(r, t) = N F ~ 1/2 J2 q m q (t)e lq r with N F being the 
number of localized spins in the ferromagnet. Consistent 
with this boson mapping, the magnetization dynamics is 
described by the following action : 14 ! 15 

S F = [ dt Y,rn-MX q {id t )]- 1 m-{t), (2) 
Jc q 

where the integration is performed along the Keldysh 
contour and the bare magnon propagator is given 



by 

X ( U \ ( X »> X «^)\ (OS 

Xqiu) ~ { 0, W ) (3) 

with the following equilibrium condition: 

X£(u) = X?( U ) = 2iImX«( W ) coth(^)) 

The retarded component of X q (uj) is given by X q (u) = 
Sq 1 ^ — uj q + iau))^ 1 where a is the Gilbert damping 
constant, and co q = jHo + to q is the magnon frequency. 
Here, 7 is the gyromagnetic ratio and ui q = D cx q 2 , where 
D cx = 2Sq J e x«s is the spin- wave stiffness constant with 
J cx and a 3 s being the exchange energy and the effective 
block spin volume. 

In the nonmagnetic metal, the dynamics of the spin 
density s can be described by the action 17 

S N = f dtY,st h {t)[Xkiid t )]- l sl{t), (5) 

where = (s? ± is?)/2 is defined by Sk = 

^w" 1/2 Ep c p+fe crc P with o"i c l = ( c P ,t' c P,;)' and Nn 
being the Pauli matrices, the electron creation operator 

for spin projection f and J,, and the number of atoms 

in the nonmagnetic metal. The equilibrium spin-density 

propagator is given by 

*W - (**$>) (0» 

with the following equilibrium condition: 

X3f(w) = [xZ(u)T, x£» = 2iIm X «(c)coth(^j)7) 

The retarded part of \ is given by^ x§ ( w ) = Xjv(1 + 
\%k 2 — iwT s f) _1 with xn, Aat, and r s f being the para- 
magnetic susceptibility, spin diffusion length, and spin 
relaxation time, the form of which is consistent with the 
corresponding diffusive Bloch equation [see Eq. (flU| be- 
low]. 

Finally, the interaction between magnons and spin den- 
sity at the interface is given by 

Sp.n = I dt S -7== m- q (t).s k (t), (8) 

JC k y/l\Fl\N 

where J^ q is the Fourier transform of J s d(r) = 
Jsd£,o(r) with J s d being the s-d exchange interaction be- 
tween conduction-electron spins and localized spins, and 

&(r) = EroGN-N interface a S S ( r ~ r o) ; 

It is instructive to point out that in the spin-wave re- 
gion and in the classical limit with negligible quantum 
fluctuations, a system described by Eqs. ©, ([5]), and © 
is equivalen t 15 ' 19 to a system described by the stochastic 
Landau-Lifshitz-Gilbert equation, 

8 t M = [ 7 (H off + h)- fya] x M + -§-M x d t M$) 
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FIG. 2: (a) System composed of ferromagnet (F) and non- 
magnetic metals (N) divided into the three temperature do- 
mains of Fi/Ni, F2/N2, and F3/N3 with their local temper- 
atures of Ti, T2, and T3. (b) Temperature profile. 



coupled with the Bloch equation^ 

d t s = (D N V 2 - t^Ss + -j^-M x a + I, (10) 



H Q z + (IW7)V 2 (M/M S ), D 



N 



where H e g 

^n/tsI is the diffusion constant, and Ss(r) = s(r) — 
s o£,o( r )M(r)/M s is the spin accumulation with the local 
equilibrium spin density so = XNSoJ s d/h~- The noise field 
h represents thermal fluctuations in F with (ft-j(r, t)) = 

and (ft i (r,t)ft J -(r',f)) = ^^^f^MC*" - r')<S(* - t')^ 
while the noise source I in N satisfies (?j(r,t)) = and 

(Z i (r,t)i j (r',t / )) = 2fcBT S xw ° 3 Mfr - r')*(t - t'W with 
the lattice constant a, both of which are postulated by 
the fluctuation-dissipation theorem. 

In this section we focus on the "local" spin injection 
from F\ into N\. The spin current induced in N\ can 
be calculated from the linear response expression of the 
magnon-mediated spin injection given in the Appendix 1X1 
[Eq. (|A4|) ]. Consider the process P\ shown in Fig. [3] (a) 
where magnons travel around the ferromagnet F\ with- 
out feeling the temperature difference between F\ and F2 . 
Using the standard rules of constructing the Feynman 
diagram in Keldysh space,— the corresponding interface 
Green's function Ck, q (oj) for the correlation between the 
magnons in F\ and the spin density in N\ [Eq. (|A4[> ] can 
be written in the form 



Xk{0j)X q (uj) 



(11) 



y/N N N F 

where Nn and N F are the number of lattice sites in Ni 
and Fx. Substituting Eq. (JTTJ) into Eq. and em- 

ploying the equilibrium conditions [Eqs. ^ and ([7))]. we 
obtain the expression for the injected spin current 



Ni 



4iV int J s 2 d 5 2 ^ 



V2fi?N N N F j£ 



1mx£{u)hnX*{u) 



coth( 



) - coth( 



2k B T, 



(12) 



where we have introduced the shorthand notation 
r _ r°° anc j ]y j s t ne number of local- 

JUJ J — OO 27T ' 1111 

ized spins at the N±-Fi interface playing a role of 
the number of channels. The oj integration can be 
performed by picking up only magnon poles under 
the condition ahLj q -C knTj\f 1 ,kBTF 1 (always satis- 
fled for YIG), giving f a Imx fe (w)ImX q (w)[coth(^)] a 

— |lmXfe (S g ) [coth( 2k ^ T )] . By making the classical ap- 
proximation coth ( y ) sa 2 ^g T , we obtain 



with the 



where Ti = da; /„* rfy p^-^^g-g— 
dimensionless variables x — fcAjv and y = Hui q /(2J cx So) 1 
and we used the relation Np 1 ^ = (27r)~ 2 J y/ydy. 



MAGNON-MEDIATED SPIN SEEBECK 
EFFECT 



III. 



Equation (JT3J) means that, through the "local" process 
Pi shown in Fig. [SJa), the spin current is not injected 
into the nonmagnetic metal iVi when F\ and N\ have 
the same temperature. That is, the "local" process can- 
not explain the experiment where no temperature dif- 
ference exists between the YIG film and the attached Pt 
film. A way to account for the experiment within the "lo- 
cal" picture is to invoke a difference between the phonon 
temperature and magnon temperature. 23 In this paper, 
on the other hand, we take a different route and consider 
the effect of temperature gradient within the YIG film on 
the spin injection into the Pt terminal. 

The basic idea of our approach is as follows. The above 
result [Eq. (fT3"|)] that the injected spin current vanishes 
when T Fl — originates from the equilibrium condi- 
tion of the magnon propagator [Eq. When magnons 
deviate from local thermal equilibrium by allowing the 
magnons to feel the temperature gradient inside the fer- 
romagnet, the magnon propagator cannot be written in 
the equilibrium form, and it generates a nontrivial contri- 
bution to the thermal spin injection. The relevant "non- 
local" process P[ is shown in Fig. [31(a) in which magnons 
feel the temperature difference between F\ and F^- The 
interaction between F\ and F% is described by the action 



S. 



F-F 



dt 



1,1' 



N F 



m q {t) ■ m_,/(t), (14) 



where J£ x q is the Fourier transform of J cx (r) = 

JUl{r) With &(r) = Er eF-F interface al 5 ^ " r 0)- 

We now regard the whole of the magnon lines appear- 
ing in the process P[ as a single magnon propagator 
8X q {bS), namely, 

5X q {u) = ^£|jr 9 ?*<zM^)* g M-(15) 
F q' 
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FIG. 3: (a) Feynman diagrams expressing the spin current 
injected from the ferromagnet (F) to the nonmagnetic metals 
(iV). The thin solid lines with arrows (bold lines without ar- 
rows) represent electron propagators (magnon propagators), 
(b) Spatial profile of the calculated spin current. 



Then the propagator is decomposed into the local- 
equilibrium part and nonequilibrium part as 2 ^ 



5X q (io) 



where 



q ~ I 0, 8X l q ^ A 



(16) 



(17) 



is the local-equilibrium propagator satisfying the local- 
equilibrium condition, i.e., SX l q ~ eq ' A = [SX l q ~ eq ' R ]* and 
8X l ~ eq > K = [8X l - eq ' R - 5X l - eq > A ] coth(3^,) with 

F q < 

while 

/n & Yn-eq,K \ 

(19) 



dx n-eg 



0, SX^~ eq - K 

0, 



is the nonequilibrium propagator with 8X q ' eq,K (uj) given 
by 



5X^ K {u) = 



\2J c l- q ' S \ 2 



X R (u)-X$(u) 



x l^«fH| 2 [coth( 5 ^-)-coth( 5 ^2)}) 

Note that the local equilibrium propagator [Eq. ([T7|] 
does not contribute to the "nonlocal" spin injection. 

When we substitute Eq. fp~6|) into Eq. (|A4j) and use 
Eq. (ITTT) with X q (uj) being replaced by SX q (uj), we ob- 
tain the following expression for the magnon-mediated 
thermal spin injection: 



jNl = Z 4J s 2 d ^o(2J cx 5 ) 2 7V i ntiV/ 



q,q',k ' 



y/2h 2 N F N N 

x |A*H| 2 ImA*H[coth(^) - coth(^%) 



where N[ nt is the number of localized spins at the Fi- 
F 2 interface, and we used = Tp i = T,- (i = 

1,2). The u> integration can be performed as be- 
fore, giving J a lm X R (uj)\X R (uj)\ 2 lmX R (oj) [coth(, hu ^ 



hUJ a 



5(u q - uj q ,)lvax R {uJ q ) [coth(^fr-) 



*-2k B T 2 >1 ™ 2aZ q 

coth( 2 ^ CJ y 2 )], which suggests that the magnon modes 
with dirferent q's do not interfere with each other. With 
the classical approximation coth( ^ q T ) w 2 ^g T , we ob- 
tain 



rjVl 



Mnt ( J 2 d s o ) Xnt s1 (a/ A jy ) 
8\/2vr 5 fi 3 a(A/as) 



-T 2 k B ST, (22) 



where 5T = Ti - T 2 , A is the size of Fi 
along the temperature gradient, and Y 2 = 

Jo' dx Jo d V Ui+x'y+v*(2SoJ m T«/h)*] which is approx- 
imated as T 2 0.1426 (T 2 « 0.337ft/25 J cx T sf ) for 
2S J ex T s {/H<l (for 2^0 J eK T sl /H > 1). 

The spin current J^ 3 injected into the right terminal 
A^ can be calculated in the same manner by consider- 
ing the process P3, which gives I^ 3 = —I^ 1 from the 
relation Ti — T 2 = — (T3 — T 2 ). The spin current I^ 2 in- 
jected into the middle terminal N% vanishes because the 
two relevant processes (P 2 and P 2 ) cancel out. Therefore, 
we obtain the spatial profile of the injected spin current 
as shown in Fig. [3fb). Note that the effect of the spa- 
tial dependence of magnetization M[T(r)] through the 
local temperature T(r) is already taken into account in 
our treatment because the temperature dependence of 
M in the magnon region is automatically described by 
the number of thermal magnons discussed in this paper. 

For an order of magnitude estimation, we com- 
pare Eq. (|22p with the experiment. 9 By using 0# « 
0.0037^26 p = 15.6xl(T 8 fmi, w = 0.1mm, X N » 7nm, 
r sf w 1 ps, a = 2 A, a s = 12.3 A, So = 16, a « 5 X 10~ 5 , 8 
Xat = 1 x 10 _6 cm 3 /g f 21 and Ai nt = 0.1 x 4mm 2 /<z|, the 
s-d exchange coupling extracted from the previous ferro- 
magnetic resonance experiment^ (J s( j ~ lOmcV) can ac- 
count for the spin Seebeck voltage Vishe/ST a; 0.1 /xV/K 
observed at room temperature. 

Finally, we comment on the issue of length scales asso- 
ciated with the SSE. In the original SSE experiment for a 
metallic ferromagnet $i the signal maintained over several 
millimeters was a big surprise because the spin diffusion 
length for that system is much shorter than a millimeter. 
Concerning the magnon-mediated SSE in an insulating 
magnet 9 which we have discussed, it is of crucial im- 
portance to recognize that the length scale relevant to 
the SSE is related to magnon density fluctuations and is 
given by longitudinal fluctuations of magnons, while the 
magnon mean free path is related to magnon dephasing 
and is given by transverse fluctuations of magnons^ It 
was shown by Mori and Kawasaki 3 ^ that these two length 
scales do not coincide with each other since they obey 
quite different dynamics, and it was demonstrated that 
in a certain situation the length scale of magnon density 
fluctuations (which is relevant to the SSE as well) is much 



5 



longer than the magnon mean free path [see Eq.(6.33) in 
Ref . where the length scale of long- wavelength magnon 
density fluctuations is infinitely long].— The notion of 
these two different length scales is the key to understand- 
ing the length scales observed in the SSE experiment in 
an insulating magnet.— 



IV. CONCLUSION 

We have developed a theory of the magnon-mediated 
spin Seebeck effect in terms of the canonical framework of 
describing transport phenomena, i.e., the linear- response 
theory, and shown that it provides us with a qualitative 
as well as quantitative understanding of the spin See- 
beck effect observed in a prototypical magnet, yttrium 
iron garnet^ Because the carriers of spin current in this 
scenario are magnons, we can obtain a bigger signal for 
a magnetic material with a lower magnon damping [see 
Eq. (|22l) where the injected spin current is inversely pro- 
portional to the Gilbert damping constant a]. An ad- 
vantage of our linear-response formulation is that it can 
be easily generalized to a situation including degrees of 
freedom other than magnons, e.g., phonons and conduc- 
tion electrons, to describe a more complicated process in 
the case of metallic 1 and semiconducting systems,— and 
a calculation taking account of the effect of nonequilib- 
rium phonons will be reported in a future publication. 34 
A numerical approach to the SSE is also developed in 
Ref. IH. We believe that the present approach stimulates 
further research on the spin Seebeck effect. 



where c p = {c p ^,c p ^) is the electron creation operator 
for spin projection f and U p — P ' is the Fourier transform 
of the impurity potential U imp Er eim P uritie S S ( r ~ r o), 
and rj so measures the strength of the spin-orbit interac- 
tion.^ 

At the ferromagnet/nonmagnetic- metal interface, the 
magnetic interaction between conduction-electron spin 
density and localized spin is described by the s-d in- 
teraction [Eq. ([5])]. The spin current induced in the 
nonmagnetic metal N\ can be calculated as the rate of 
change of the spin accumulation in Ni, i.e., Ig(t) = 
J2reN 1 ( d tS z (r,t)} = (d t s z ko (t)) ko ^ , where (• • • ) means 
the statistical average at a given time t, and s& = 
^/NpjSk with s being denned below Eq. ([5]). 

The Heisenberg equation of motion for s z k a gives 



E 

q.k 



^2N F N N h\ 



-ql S k ! S k \ 



m_ q [s k ,s ho 



q.k 



y/2N F N N h 



±2s fc+fc" and 



where we have used the relation [s| , s^, 
neglected a small correction term arising from the spin 
orbit interaction assuming that the spin-orbit interaction 
is weak enough at the neighborhoods of the interface. 
Then, the statistical average of the above quantity gives 
the following spin current: 
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Appendix A: Linear- response expression of 
magnon-induced spin injection 

The Gaussian action for conduction electrons in the 
nonmagnetic metal JV^ (i — 1, 2, 3) is given by 

S N = J dt y^4(t)|u? t - (e P 5 P ,p' 
c p,p' 

+ U p - p/ [1 + lT] so (T • (p x p')]) jcp< (t), (Al) 



7j (*) 



(A3) 



where C^ q (t,t') = — i(ml q (t')s^"(t)) is the interface 
Green's function. In the steady state, the Green's func- 
tion C k (t, t') depends only on the time difference t — t' 



as C< q (t f) = /. 
the representation^ C = 
U C K -C R + C A ], we finally obtain 



&C< fc ( W )e 

' c K 

C A 



-iw(t-t') < Adopting 
and using C < = 



7^ = 



k \ 2.V/ A v// 



^ReC^H (A4) 



for the spin current I^ 1 in a steady state. As in the case 
of tunneling charge current driven by a voltage differ- 
ence^ the spin current I^ 1 can be calculated systemat- 
ically. 
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